Quantum Nonlinear Sigma Model for Arbitrary Spin Heisenberg Antiferromagnets 
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In this Letter, we derive a quantum nonlinear sigma model (QNLSM) for quantum Heisenberg 
antiferromagnets (QHA) with arbitrary S (spin) values. A upper limit of the low temperature is 
naturally carried out for the reliability of the QNLSM. The S dependence of the effective coupling 
constant and the spin wave velocity in the QNLSM are also obtained explicitly. The resulting 
spin wave velocity for 2-dim spin- 1/2 QHA highly concurs with the experimental data of high T c 
compound La2CuC>4. The predicted correlation lengths for 2d QHA and spin-gap magnitudes for 
Id QHA also agrees with the accurate numerical results. 
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Quantum magnetism in low dimensional strongly cor- 
related systems is a central issue in modern condensed 
matter physics. For example, the parent compounds 
of the cuprate high T c superconductors are essentially 
antifcrromagnetic Mott insulators described by a two- 
dimensional spin-i QHA pj. Such a low-dimensional 
QHA system has not been fully solved quantum mechan- 
ically. Only in the large S limit has Haldane shown that 
the lattice QHA can be described by the QNLSM 0. 
Since then, the QNLSM has become a good candidate for 
the phenomenological description of the low-dimensional 
QHA at low temperatures for various S values 0, 0, 0. 

However, two crucial questions remain for the QNLSM 
approach: (1) Why are the predictions of QNLSM ob- 
tained in the large-.!? limit coincidentally consistent with 
experimental data of the low energy QHA with small 
S values? and (2) What is the upper limit of the low 
temperature for the reliability of QNLSM? These two 
questions, along with other difficulties for the QNLSM, 
have been addressed in some recent literatures 6J. But 
they have not yet been satisfactorily solved in a simple, 
consistent approach. In this Letter, with the topologi- 
cally invariant spin variable path integral approach, we 
resolve these problems by deriving a QNLSM from the 
lattice QHA for arbitrary S values. 

Let us first briefly recall the large-S" approach of the 
QNLSM. The partition function of a spin system is usu- 
ally expressed in terms of the spin variable path integral 
as follows (in the unit h= Ub = 1): 



Z = Tre 



-0H 



J T>[n]eJf 



dT{iSA-n-(n\H\n)} 



(1) 



where the first term in the exponent iS J drA. ■ fl = 

13 dr{fl\4-\n) = iSu>(fl) is a topological Berry phase 
and A is a U(l) monopole potential; |C1) is a spin 
coherent state while $7 is a unit vector along which the 
spin operator with spin quantum number 5* is maximally 
aligned in |fl). For the lattice Heisenberg model (HM). 



H = JJ2<ij> S i ■ S j (J>0), Eq.(|TJ becomes 
Z H = f^n^dr^SE^^-js^E^^-^}^ (2) 



By minimizing H(Q) = JS 2 ^2 <i j > Qi ■ flj, one can 
find the classical ground state (Neel state) which sponta- 
neously breaks the SO (3) symmetry. Then by expanding 
the action around the ground state, the spin- wave theory 
(SWT) of HM, which describes the long wavelength spin 
modes ;8j, can be easily derived. 

But, according to the Mermin- Wagner's theorem^, no 
symmetry can be spontaneously broken in one- or two- 
dimensional HM for a finite T(> 0). To derive an effec- 
tive long wavelength action that retains the full spin ro- 
tational symmetry, Haldane considered the large-S* limit. 
In the large-S 1 limit, the path integrals of Eq.J2J) are 
dominated by the semiclassical equation: iSfl x ft = 
■ By separating the semiclassical solution tti into 
a slowly varying Neel order unit vector (— l)'n(xj) plus 
a slowly varying magnetization density field perpendic- 
ular to n(xi) (Haldane's mapping), and then taking the 
continuous limit and integrating out the magnetic density 
field, Haldane shows that Eq.J2J is reduced to a QNLSM, 



Z H = / V[n] 



J2irS0in] -(A d - 1 /2g ) f d d+1 x d^tid^l 



(3) 



defined in the <i+l-dimensional space (ar, 
x d+l) = (a;i ? ... ,x d ,c T). Where, g a = 2y/d/S is a di- 
mensionless coupling constant, c D = 2y/dJSa the spin 



wave velocity and A 



the inverse of the lattice 



spacing. The imaginary time (temperature) variable r 
ranges from to j3 = 1/T. The exponent 0[n] in Eq.© 
is a topological factor associated with the Berry phase. 

However, as we will elaborate, the derivation of the 
QNLSM based on the large-S* expansion should be im- 
proved from the very beginning. As a long-standing 
problem in the construction of generalized phase space 
path integrals^^i the Eq.Q is not well defined. The 
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main problem arises from the assumption, used in deriv- 
ing Eq.Q, that |Cl(r + St)) - |«(t)) is of order 0(St). 
Although this assumption has been widely used in the 
application of generalized phase space path integrals, it 
has never been justified[lO|. In fact, in effective field the- 
ories, there always exist simultaneous rapidly and slowly 
varying paths in the path integral formalism that are as- 
sociated with short and long range quantum fluctuations, 
respectively. The effective action for slowly varying mo- 
tions can be properly obtained by integrating over short 
range quantum fluctuations |l 1| . However, in Eq.Q), only 
slowly varying motions are retained; the short range 
quantum fluctuations have been simply ignored. 

To overcome the shortcomings involved in the deriva- 
tion of Eq.QJ, we begin with the discrete form of the par- 
tition function obtained exactly from the coherent state 
representation [T^ : 



N 



N 



Z= lim TT ^(O fc )cxp{V In(r2 fc |f2 fc - 1 

N^oo - 1 -*- z — ' 
k=l k=l 

(n k \H\n k - 1 



(ft k \ft 



fc-i\ 



}, (4) 



where \fl N ) = \ft°) because of the periodicity of the 
trace), e — (3/N is infinitesimal as iV — * oo. The 
slowly varying motion means that |Oj) — |fij_i) varies 
smoothly in the interval e such that it can be written 
as a time derivative e\ft). The rapidly varying motion of 
|Oj) — in the interval e is related to the short range 

quantum fluctuations, we label it as 6 ft. The assumption 
of \ ft k ) — |fJ fe_1 ) being of order 0(e) only keeps the slowly 
varying motions ft, while the short range fluctuations 6ft 
is ignored. Thus, the continuous-time limit of Q results 
in the conventional spin variable path integral of Eq. • 
To include the contribution of the short range fluc- 
tuations, one should expand the nearby coherent state 
overlap to the second order terms that either are ex- 
clusively slowly varying motions or include at least one 
rapidly varying motion. The topologically invariant 

terms limAr-^oo Ylk=i ^ n (^k\ftk-i) can then be uniquely 
expressed by 







SJ eZrjiA • ft + ifl ■ (ft X Sfl) 6ft ■ Sfl}, (5) 

where the parameter t\ = 1/T\ is an intrinsic short- 
est time scale (an upper limit of low temperature) for 
distinguishing the slowly varying and rapidly varying 
motions. We will later discuss this timescale in de- 
tail. The second and third terms in Eq.10 are usually 
ignored in the conventional derivation of phase space 
path integrals. For the Hamiltonian term in Eq.(0}, 

limjv^oo J2k=i 6 (n^n 1 *- 1 ) > since lt is alread y propor- 
tional to e, we only keep the off-diagonal expansion up 
to the quadratic order of 6ft: 



(It 



[H[fl] 



dH[ft] 
dft 



6ft 



d 2 H[ft] 



dfl n dft, 



-6fl a 5fl a ,y (6) 



where a, a' are indices of spin components. Substituting 
Ecis. I|5I6[) into Eq.Q}, one gets 

Z = J V[ft\D[6ft] exp J dT^iSA-fl- H[ft] 

dH[ft] 



iSfl x ft 



dft 

l Saa , + d2H ^ 
ta aa dfl a dfl a > 



■ 6fl 

6fl a 5fl a ,y (7) 



which describes both the slowly varying motion ft and 
the short range fluctuations 6ft. 

If one takes the semiclassical limit, iSft x ft — 
0. Eq.Q is simply a variation expansion of the path inte- 
gral with respect to the semiclassical dynamics Hal- 
dane used[2], except for the geometrical term ^- which 
cannot appear in Haldane's mapping. However, we must 
emphasize that Eq. Q is derived by carefully treating the 
off-diagonal elements of nearby coherent states in Eq. Q 
in terms of the short range fluctuation 5 ft and the slowly 
varying motions ft. Since there is no semiclassical expan- 
sion to begin with, it is not necessary to take the semi- 
classical limit by letting the second term vanish. Instead, 
one can integrate out the short range fluctuations 6 ft and 
obtain a low energy effective action for the long wave- 
length spin modes. Since no semiclassical approximation 
is made in this procedure, the resulting long wavelength 
effective action should be valid for arbitrary S values. 

Next, we apply Eq.© to the HM. To specify the an- 
tiferromagnetic ordering, let the slowly varying fti = 
(— l) l n(xj), here the Need order n(xj) is a unit vector 
|n(xj)| = 1. Then, by taking the space continuous limit 
J2i —>■ -^s J d d x where a is the lattice spacing: 

H[ft] = js 2 n< • 

<«?> 

JS 2 



-dJS 2 N 



2a d ~ 2 



d d xY,[dMx) ■ d k n(x)}, (8) 



fe=i 



dH[ft] 
dft 
d 2 H[fl] 



■6ft^0, 



(9) 



5ft a 5ftf3 — — I d d x6fl(x) ■ 6fl{x). (10) 



dfl a dflp 

substituting (|8I10() into (J7J) , and integrating out the short 
range fluctuation 6 ft, we obtain, 



Z H = f V[n] e ^se[n] exp | 



,1-d 



2.9, 



d(c s T) 
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This is a QNLSM of the low energy QHA for arbitrary 
S values, where the coupling constant and spin-wave ve- 
locity are given by 



2SJ 



2JSa\ d 



2SJ 



(12) 
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The topological phase factor 27r50[n] = 
2irS^2 i (~l) l A(xi) ■ h(xi) remains the same as in 
Haldane's derivation. 

In the large- 5 limit for fixed T\, g s — * 2y/d/S,c s — > 
2y[dJSa. This reproduces the large-5 QNLSM. The 
difference between (J3J) and primarily comes from 
the contribution of the short range fluctuations, the 
1 /ta = T\ term in (JJJ which cannot be included in Hal- 
dane's mapping 0. However, this term plays an impor- 
tant role in the derivation of a consistent semiclassical 
dynamics [l3j. Indeed, T\ is an upper limit of the low 
temperature scale for the reliability of the QNLSM. Usu- 
ally one assumes that there should be no intrinsic cut- 
off for the imaginary time variable r because quantum 
fluctuations exist on all time scales in path integrals Q. 
But a low energy effective theory constructed from path 
integral is defined by integrating over high energy quan- 
tum fluctuations above certain energy scale 01 • Without 
such an intrinsic cutoff, namely, let ta — * (Ta — > oo), 
Eci. (fTT|l reduces to 



Z H oc / D[n] e 



i27rSe[nl 



exp 



r _ ps_ 

I 2T 



d d x |V x n| : 



]}• 

(13) 

where p s — JS 2 a 2 d is the spin stiffness. Except for 
the topological phase, this is the classical d-dimensional 
NLSM rather than a quantum d + 1-dimensional NLSM 
that Haldane obtained^]. This is because in the limit 
ta — > 0, the strong canonical fluctuation in Eq.© smears 
the dynamical fluctuation of Ea. (|l(Jfl so that only the 
classical Hamiltonian Eq.© remains. 

Now, let us discuss how to consistently determine this 
timescale. The lattice spacing a indicates the existence 
of an intrinsic momentum cutoff A in the ci-dimensional 
momentum space: A = 2^/n[T(d/2 + l)] 1/d /a = L/a. 
Using the energy-momentum relation of the spin wave, 
E = c s k, one can find the intrinsic energy cutoff (the 
inverse of the shortest time scale ta) Ta = c s A/27r[j|. 
Combined with i|12fl . we get 



Ta 
J 



SL 2 
Air 2 



I6ir 2 d 



(14) 



For d = 2 and 5 = 1/2, we have L = 2^/tt and thus 
T\/ J ~ 0.97. This determines quantitatively a low tem- 
perature upper limit for the reliability of QNLSM to the 
2-dim spin-1/2 QHA: < T/J < T A /J ~ 1.0. Mean- 
while, the spin wave velocity c s can also be explicitly de- 
termined from Ea. l|12|l and Eq. i|14fl . For La 2 Cu04, which 
is a typical 2-dim spin-1/2 QHA with a — 3.79A and J ~ 

1500if, we obtain (keeping K) hc s — 2JSa^J d + j^y ~ 

0.85 eV A. This is in excellent agreement with the exper- 
imental data hc s = 0.85 ± 0.03 eV A0- 

Our main results, i.e., Eas. (|12|) and 114|) . can be further- 
tested against the known results for the 2d QHA, among 
which the quantum Monte Carlo date is almost exact. 
Up to the three-loop correction, the QNLSM predicts 3 
the asymptotic scaling behavior of the correlation length 



in the renormalized classical regime as 

C 3i =^exp(l/<)[l-O.5i + 0(i 2 )], 



(15) 



where, A — | is a temperature-independent pre- 
f actor, t — 2^jf- is the dimensionless temperature, and c s 
and p s are the renormalized spin-wave velocity and spin 
stiffness which can be consistently determined by large-5 
expansion. The predicted formula £3/ is extremely sensi- 
tive to the spin stiffness and is consistent with the QMC 
at very large correlation lengths (low temper- 
atures) for 5—1/2 when the best-fit value p s — 0.1800 is 
used. However, at moderate correlation lengths, highly 
accurate QMC data0| and series expansionsjl^ indicate 
a significant discrepancy which rapidly increases with 5. 
The asymptotic scaling at the three-loop sets in at cor- 
relation lengths larger than 10 5 for 5 = 1/2^^ and cos- 
mological lengths for larger 5. 

Quite strikingly, according to the basic assumption of 
the large-5 approach, the discrepancy between the the- 
ory and numerics, if it exists, should be significant only 
for small 5. By contrast, in our derivation, the asymp- 
totic scaling behavior Eq. l|15|l holds for arbitrary 5 at 
low temperatures, provided the effects of the intrinsic 
scale Tv are correctly taken into account. Note that the 
temperature dependence in Eq. I|15|) comes from a simple 
assumption that at nonzero temperature T the correla- 
tion length is much larger than the finite extent c s /T 
along the Euclidean time direction in the renormalized 
classical regime [3J, In our case, the time extent is re- 
placed by c s (l/T— 1/Tv), making the assumption more 
reasonable. Therefore, we only need to replace T by the 
re-scaled temperature T = if^f in Eq.JTHJl. The valid- 
ity of this simple re-scaling requires that T < 2np s , or 
T < 0.5 J for 5 = 1/2. Interestingly, the re-scaling does 
not change the two-loop asymptotic scaling behavior. On 
the other hand, the shifts in c s and g s only modify the 
pre- factor A but keep p s unchanged. Figure ^ shows 
the deviations of various results from two-loop asymp- 
totic scaling exp(l/t) as a function of t. The three-loop 
results of Ea. (|15fl in terms of t and t — are plot- 

ted as the 3-loop old (dashed) line and the 3-loop new 
(solid)line, respectively. For our three-loop result, the 
scaling regime begins at roughly £ sw 10 2 — 10 3 , when 
T ps 0.3J. As a comparison, a suggested four-loop (dot- 
ted) line^jjj is plotted by adding the correction 0(t 2 ) in 
Eci. (|15|) before re-scaling. The coefficient of this term is 
—0.75 by fitting the highly accurate QMC data which 
is apparently too lar ge t o be obtained within reasonable 
four-loop corrections [15j. 

Our results can also be tested against the known re- 
sults for the Id QHA. Note that in Id, the topological 
term plays a crucial role, leading to the quantum crit- 
ical and disordered phases for half integer and integer 
spins, respectively 0, 0. In the disordered phase, T\ 
can be determined by using the relation E 2 = c 2 k 2 + A 2 . 
This only leads to a correction C((^j) 2 ) to Eq.JTl}. 
For 5 = 1/2,1,3/2 and 2, the accurate numerical re- 
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FIG. 1: Deviations of various results from two-loop asymp- 
totic scaling as a function of t. The QMC result is obtained 
from Ref.ll5ll. 



suits for the values of c s /(aJ) obtained by the Density 
Matrix Renormalization GroupfDMROfTa. I20I l2l|| are 
1.57,2.49,3.87 and 4.65 respectively, which show a sys- 
tematic deviation from 25. While, by Ea. l|12[l . they are 
1.28, 2.55, 3.84 and 5.00 respectively, providing much bet- 
ter predictions. The Id QNLSM also predicts gap[i"s| in 
the disordered phase as A s = Bc s exp(— 2ir/g s ), where 
B is a ^-independent fitting parameter of order of 1. 
The DMRG results are A, = 0.411J for 5" = l|!2; 
A 2 = 0.085 J for S = 2^ |23|. By using the best- 
fitting parameter, B ss 2, we find that our derivation 
gives Ai = 0.438 J and A 2 = 0.076 J, respectively, while 
the large-5* approach gives Ai = 0.172 J and A 2 = 0.015, 



respectively. Therefore, the present Id QNLSM provides 
a better approximation for the spin-gap magnitudes in 
the Id integer QHA. 

In conclusion, by using a topologically invariant spin 
variable path integral approach, we resolve the two cru- 
cial problems in the QNLSM description of QHA as men- 
tioned in the beginning of this Letter. The basic param- 
eters in the QNLSM are unambiguously defined for ar- 
bitrary spin values. The primary tests discussed above 
show that the quantum fluctuations in the QHA, which 
are usually underestimated in the large-.!? approach, are 
now more properly described. It should be emphasized 
that the construction of a low energy effective field theory 
from the extended phase space path integrals developed 
in this Letter is a general approach, in which the shortest 
time scale plays a crucial role for self-consistency. This 
approach can be applied to other generalized phase space 
path integrals [12] for the study of low energy physics in 
various strongly correlated systems. 
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